Vacuum-field Rabi oscillations in atomically doped carbon nanotubes 
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We report a strictly non-exponential spontaneous decay dynamics of an excited two-level atom 
placed inside or at different distances outside a carbon nanotube. This is the result of strong non- 
Markovian memory effects arising from the rapid frequency variation of the photonic density of 
states near the nanotube. The system exhibits vacuum-field Rabi oscillations when the atom is 
close enough to the nanotube surface and the atomic transition frequency is in the vicinity of the 
resonance of the photonic density of states. 
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It has long been recognized that the spontaneous emis- 
sion rate of an excited atom is not an immutable prop- 
erty, but that it can be modified by the atomic envi- 
ronment. Generally called the Purcell effect the phe- 
nomenon is qualitatively explained by the fact that the 
local environment modifies the strength and distribution 
of the vacuum electromagnetic modes with which the 
atom can interact, resulting indirectly in the alteration 
of atomic spontaneous emission properties. 

The Purcell effect took on special significance recently 
in view of rapid progress in physics of nanostructures. 
Here, the control of spontaneous emission has been pre- 
dicted to have a lot of useful applications, ranging from 
the improvement of existing devices (lasers, light emit- 
ting diodes) to such nontrivial functions as the emis- 
sion of nonclassical states of light P|- In particular, the 
enhancement of the spontaneous emission rate can be 
the first step towards the realization of a thresholdless 
laser 3] or a single photon source Q. The possibility to 
control atomic spontaneous emission was shown theoret- 
ically for microcavities and microspheres 0, optical 
fibers photonic crystals semiconductor quantum 
dots (2] ■ Recent technological progress in the fabrication 
of low-dimensional nanostructures has enabled the exper- 
imental investigation of spontaneous emission for micro- 
cavities photonic crystals quantum dots [l2] . 

Recent successful experiments on encapsulation of sin- 
gle atoms into single-wall carbon nanotubes (CNs) M| 
and their intercalation into single-wall CN bundles [iJI 
stimulate an analysis of atomic spontaneous emission 
properties of such systems. Typically, there may be two 
qualitatively different regimes of the interaction of the ex- 
cited atomic state with a vacuum electromagnetic field in 
the vicinity of the CN. They are the weak coupling regime 
and the strong coupling regime. The former is charac- 
terized by the monotonous exponential decay dynamics 
of the upper atomic state with the decay rate altered 
compared with the free-space value. The latter is, in con- 
trast, strictly non-exponential and is characterized by re- 
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versible Rabi oscillations where the energy of the initially 
excited atom is periodically exchanged between the atom 
and the field. In this Letter, we develop the quantum the- 
ory of the spontaneous decay of a two-level atom near a 
CN, derive the evolution equation of the upper state of 
the system and, by solving it numerically, demonstrate 
the strictly non- exponential spontaneous decay dynamics 
in the case where the atom is close enough to the CN sur- 
face. In certain cases, the system exhibits vacuum-field 
Rabi oscillations - a result already detected for quasi- 
2D excitonic and intersubband electronic transitions in 
semiconductor quantum microcavities [isl Il6| and never 
reported for atomically doped CNs so far. 

The quantum theory of the spontaneous decay of ex- 
cited atomic states in the presence of a CN involves an 
electromagnetic field quantization procedure. Such a pro- 
cedure faces difficulties similar to those in quantum optics 
of 3D Kramers-Kronig dielectric media where the canon- 
ical quantization scheme commonly used does not work 
since, because of absorption, operator Maxwell equations 
become non-Hermitian As a result, their solutions 
cannot be expanded in power orthogonal modes and the 
concept of modes itself becomes more subtle. We, there- 
fore, use an alternative quantization scheme developed 
in Ref. 0|, where Fourier- images of electric and mag- 
netic fields are considered as quantum mechanical ob- 
servables of corresponding electric and magnetic field op- 
erators. The latter ones satisfy the Fourier-domain op- 
erator Maxwell equations modified by the presence of 
a so-called operator noise current written in terms of a 
vector bosonic field operator and an imaginary dielec- 
tric permittivity. This operator is responsible for cor- 
rect commutation relations of the electric and magnetic 
field operators in the presence of medium-induced ab- 
sorbtion. The electric and magnetic field operators are 
then expressed in terms of a continuum set of the bosonic 
fields by means of the convolution of the operator noise 
current with a classical electromagnetic field Green ten- 
sor of the system. The bosonic field operators create and 
annihilate single-quantum electromagnetic medium exci- 
tations. They are defined by their commutation relations 
and play the role of the fundamental dynamical variables 
in terms of which the Hamiltonian of the composed sys- 
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tern " electromagnetic field + dissipative medium" is writ- 
ten in a standard secondly quantized form. 

We consider a two-level atom positioned at the point 
ta near an infinitely long single-wall CN. The atom in- 
teracts with a quantum electromagnetic field via an elec- 
tric dipole transition of frequency loa- The atomic dipole 
moment, d = d^e^, is assumed to be directed along the 
CN axis assigned by the unit vector of the orthonor- 
mal cylindric basis {6^,6^,6^}. The contribution of the 
transverse dipole moment orientation is suppressed be- 
cause of the stron g de polarization of the transverse field 
in an isolated CN [13 . The Hamiltonian of the system is 
given by 

H = J dnj^ duj fiuj p (R, tj)/(R, ^) + ^ Tt-oja cfz 

-[at^W(rA)d.+h.c.] (1) 

with the three terms representing the electromagnetic 
field (modified by the presence of the CN), the two- 
level atom and their interaction (within the framework 
of electric dipole and rotating wave approximations |^), 
respectively. The operators /^(R, lj) and /(R, w) are the 
scalar bosonic field operators defined on the CN surface 
(R = {Ren, (f'l Z} is the radius-vector of an arbitrary 
point of the CN surface). They play the role of the fun- 
damental dynamical variables of the field subsystem and 
satisfy the standard bosonic commutation relations. The 
Pauli operators az, <J, ct^ describe the atomic subsystem. 
The operators E(='=^(r^) represent the electric field the 
atom interacts with. For an arbitrary r — {r, ip, z}, they 
are defined as E(r) E(+)(r) + E(-)(r) with 

M+\r)= F,{r,io)dLu, E(-)(r) = [EW(r)]^ (2) 
Jo 

Here, E satisfies the Fourier-domain Maxwell equations 
V X E(r,tj) = ikH{r,uj), 

V X H(r, Lu) = -ifcE(r, uj) + — i(r, uj), (3) 

c 

where H stands for the magnetic field operator [defined 
by analogy with Eq. (0)], k — uj/c, and 

i(r, uj) = y"dR (5(r-R) J (R, u;) = 2J(i?e„ ,^,z, w)(5(r-i?,„) 

(4) 

with J = ^yhuJKeazz (R, w) / 7r/(R, 10)6^ being the oper- 
ator noise current density associated with CN-induced 
absorbtion {azz is the CN surface axial conductivity per 
unit length; in describing CN electronic properties we use 
the model of the axial conductivity ^20\). 
From Eq. Q it follows that 

E(r,a;) fc /dRG(r,R,w)-J(R,w) (5) 



[and H = {ik)~^V x E accordingly], where G(r, R, uj) is 
the Green tensor of the classical electromagnetic field in 
the vicinity of the CN. Its components satisfy the equa- 
tion 

J2 (VxVx-fc2)^^G„,(r,R,c^) = 5(r-R), (6) 



together with the radiation conditions at infinity and the 
boundary conditions on the CN surface. The latter ones 
are obtained from the classical electromagnetic field bou- 
ndary conditions derived in Ref. 20] by means of sim- 
ple relations valid for r^^r^: Ea{r,uj) ^ ikGaz{j^,J^A,^) 
and Ha = {'ik)~^J2f3 7 ^aPi'^ P^i {^aPi is the totally anti- 
symmetric unit tensor). The Hamiltonian along with 
Eqs. ©-(El is the modification of the Jaynes-Cummings 
model 17] for an atom in the vicinity of a solitary CN. 

When the atom is initially in the upper state and the 
field subsystem is in vacuum, the time-dependent wave 
function of the whole system can be written as 



lV(t)) =C„(t)e- 



-j(wA/2)t 



u)m) 



+ dv d^G(r,c.,i)e-(— -/2)*lZ)l{l(r,a;)}), (7) 



where \u) and \l) are the upper and lower atomic states, 
respectively, [{0}) is the vacuum state of the field subsys- 
tem, ]{l(r,Li;)}) is its excited state where the field is in a 
single-quantum Fock state. In view of Eqs. ||2J) and lO, 
the Schrodinger equation with Hamiltonian |^ and wave 
function ||7J| yields the following evolution law for the oc- 
cupation probability amplitude C« of the upper state 



Cuir) = 1 + 



rK{r 
Jo 



^Cu{T')dT' , 



(8) 



Kir r') = ^^4^ / d..^a-)'— 

(9) 

Here, x = fiuj/2'-fo and r — 2'-fot/h are the dimensionless 
frequency and time, respectively, with 70 = 2.7 eV being 
the carbon nearest neighbor hopping integral appearing 
in the CN surface axial conductivity in the noise current 
in Eq. |0J, ^{x) = T{x)/ro{x) is the relative density (with 
respect to vacuum) of photonic states near the CN with 
r(a;) = {8'Kdl/hc'^){2jox/h)HmGzz{'rA,'rA,x) being the 
rate of the exponential spontaneous decay of the atom 
near the CN and ro(a;), its vacuum counterpart, given by 
the same formula with Im Gl^{x) = {l/6Trc){2jQx/h). The 
function ^(x) was first derived and analyzed in Ref. |l8j ]. 
Eq. (jSJ is a well-known Volterra integral equation of the 
second kind. In our case, it describes the spontaneous de- 
cay dynamics of the excited two-level atom in the vicinity 
of the CN. All the CN parameters that are relevant for 
the spontaneous decay are contained in the classical elec- 
tromagnetic field Green tensor in Eq. 0. 

In the case where the Markovian approximation is ap- 
plicable, or, in other words, when the atom-field coupling 
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FIG. 1: (Color online) Relative density of photonic states (a) 
and upper-level spontaneous decay dynamics (b) for the atom 
in the center of different CNs. 



strength is weak enough for atomic motion memory ef- 
fects to be insignificant, so that they may be neglected, 
the factor {exp[— — xa){t — r')] — l}/i{x — xa) in 
Eq. may be replaced by its long-time approxima- 
tion —Trd{x — Xa) + iV I (x — xa) {V denotes a principal 
value). Then, K{t - t') = ~hT{xA)l'^lo + i^{xA) with 
A(a;A) = [?iro(a;A)/47rx^7o] V dx x'^£,{x) / [x-xa), and 
Eq. lO yields (7u(t) =exp[(-^r/47o-l- jA)r] — the expo- 
nential decay dynamics of the [shifted by A(xa)] upper 
atomic level with the rate T{xa)- This case was analyzed 
in Ref. 

To take non-Mar kovian effects into account, we have 
solved Eqs. ^ and Q numerically. The exact time evo- 
lution of the occupation probability |(7„(t)|^ of the up- 
per state was obtained for the atom placed in the cen- 
ter inside and at different distances outside achiral CNs 
of different radii. The factor ^(x) in Eq. Q was com- 
puted in the same manner as it was done in Ref. |l8j |. 
The free-space spontaneous decay rate was approximated 
as Tq(xa) ~ (X^'^lo XA I Ti with a = 1/137 being the fine- 
structure constant plj . 

Figure^^a) presents £,{x) for the atom in the center of 
the (5,5), (10,10) and (23,0) CNs. It is seen to decrease 
with increasing the CN radius, representing the decrease 
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FIG. 2: (Color online) Relative density of photonic states (a) 
and upper-level spontaneous decay dynamics (b, c) for the 
atom at different distances outside (9,0) CN. 



of the atom-field coupling strength as the atom moves 
away from the CN surface To calculate \Cu{t)\^ in 
this particular case, we have fixed xa = 0.45 (indicated by 
the vertical dashed line), firstly, because this transition 
is located within the visible light range 0.305 < a; < 0.574, 
secondly, because this is the approximate peak position of 
^{x) for all the three CNs. The functions |C„(t)P calcu- 
lated are shown in Figure ^b) in comparison with those 
obtained in the Markovian approximation yielding the 
exponential decay. The actual spontaneous decay dynam- 
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ics is clearly seen to be non-exponential. For the small 
radius (5,5) CN, Rabi oscillations are observed, indicat- 
ing a strong atom-field coupling regime related to strong 
non-Markovian memory effects. With increasing the CN 
radius, the decay dynamics approaches the exponential 
one with the decay rate enhanced by several orders of 
magnitude compared with that in free space. 

Figure |2Ia) shows ^(x) for the atom outside the (9,0) 
CN at the different distances from its surface. The verti- 
cal dashed lines indicate the atomic transitions for which 
the functions \Cu{t)\'^ in Figures |2Jb) and|2Ic) were cal- 
culated. Both transitions belong to the visible light range, 
however xa = 0.33 is the position of a peak (for the 
shortest atom-surface distance) while xa = 0.52 is the 
position of a dip of £,{x). Very clear underdamped Rabi 
oscillations are seen for the shortest atom-surface dis- 
tance at a;^ = 0.33 [Figure [^fb)], indicating strong atom- 
field coupling with strong non-Markovity. For xa = 0.52 
[Figure [2tc)], though ^(0.52) is comparatively small, the 
spontaneous decay dynamics is still non-exponential, ap- 
proaching the exponential one only when the atom is far 
enough from the CN surface. Similar to what takes place 
in photonic crystals Q], this is because of the fact that, 
due to the rapid variation of S,{x) in the neighborhood 
of this frequency, the correlation time of the electromag- 
netic vacuum is not negligible on the time scale of the 
evolution of the atomic system, so that atomic motion 
memory effects are important and the Markovian approx- 
imation in Eq. Q is inapplicable. 

The effects predicted will yield an additional structure 
in optical absorbance/reflectance spectra (see, e.g., 0) 
of atomically doped CNs in the vicinity of the energy of 
an atomic transition. Weak non-Markovity of the spon- 
taneous decay (non-exponential decay with no Rabi os- 
cillations) will cause an asymmetry of an optical spec- 
tral line-shape (similar to exciton optical absorbtion line- 
shape in quantum dots [l^l). Strong non-Markovity of 



the spontaneous decay (non-exponential decay with fast 
Rabi oscillations) originates from strong atom-vacuum- 
field coupling with the upper state of the system splitted 
into two "dressed" states. This will yield a two-component 
structure of optical spectra similar to that observed for 
excitonic and intersubband electronic transitions in semi- 
conductor quantum microcavities . 

To conclude, we predict the existence of freque ncy 
ranges where, similar to semiconductor microcavities |l6| 
and photonic band-gap materials 0, CNs qualitatively 
change the character of atom-electromagnetic-field inter- 
action, yielding strong atom-field coupling — an impor- 
tant phenomenon necessary, e.g., for quantum informa- 
tion processing The present Letter dealt with the 
simplest manifestation of this general phenomenon — 
vacuum-field Rabi oscillations in the atomic spontaneous 
decay dynamics near a single-wall carbon nanotube. Such 
a strictly non-exponential decay dynamics gives place to 
the exponential one if the atom moves away from the 
CN surface. Thus, the atom-field coupling strength and 
the character of the spontaneous decay dynamics, respec- 
tively, may be controlled by changing the distance be- 
tween the atom and the CN surface by means of a proper 
preparation of atomically doped CNs. This opens routes 
for new challenging nanophotonics applications of atom- 
ically doped CN systems as various sources of coherent 
light emitted by dopant atoms. We emphasize that simi- 
lar manifestations of strong atom-field coupling may oc- 
cur in many other atom-electromagnetic-field interaction 
processes in the presence of CNs, such as, e.g., dipole- 
dipole interaction between atoms by means of a vacuum 
photon exchange |2^ , or cascade spontaneous transitions 
in three-level atomic systems p5| . 
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